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Abstract 

We obtain the radius of convergence of the small-amplitude approximation to the 
period of the nonlinear oscillator x + (1 + x 2 )x = with the initial conditions 
x(0) = A and x(0) = and show that the inverted perturbation series appears to 
converge smoothly from below. 

There has recently been great interest in the study of the period of the non- 
linear oscillator 

x{t) + [1 + x(t) 2 ]x(t) = 

x(0) = A, x(0) = (1) 

as a function of the amplitude A. Apparently, it aroused from the fact that the 
first-order harmonic balance method yielded the approximate frequency [1] 

= , 2 (2) 
that is not denned for A > 2. 
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By straightforward analysis of the dynamical trajectories in the x — y plane, 
where y = x, Beatty and Mickens [2] concluded that such a restriction is 
merely an artifact of the harmonic balance method. 



Later, Mickens [3] derived an explicit expression for the period 

where u = x/A. By means of this expression he proved that dT/dA < and 
obtained upper and lower bounds to the period. 

Kalmar-Nagy and Erneux [I] derived the behaviour of the period for small 
and large values of A 



T(A)^, A»l (4) 

as well as most interesting approximations to the periodic orbits in both limits. 
In particular, they showed that the trajectory u(t) satisfies the equation 



dV 
u + — = 0, 
du 

V(u) = , p = A 2 (5) 

that leads to the same expression for the period ([3]) derived earlier by Mickens 



The results of those authors clearly show that the period T(A) does not ex- 
hibit singular points for real values of A but they do not explain why the 
harmonic balance fails as shown in equation © P]. A possible explanation 
is that the harmonic balance is reflecting a singular point in the complex bi- 
plane. If it exists, then the small-amplitude expansion will have a finite radius 



2 



of convergence. 



In order to derive the small-amplitude expansion we change the integration 
variable in equation (E]) to u = cos 9 so that the period becomes 



tt/2 

T(P) = 4 / - 



dO 



F(psin 



(6) 



where 



F(z) 



e z - 1 



h u + 



(7) 



If we substitute the expansion 
1 



c i ; 



F(z) j=0 



z z 2 z 3 

1 1 1 h 

4 96 384 



into the equation (jSj) we obtain as many coefficients as desired of the small- 
amplitude series 



T(A) = T + T lP + T 2 p 2 + ... 

= Z7T 1 h • • 

\ 8 256 



(9) 



The function l/\/e z — 1 has two complex-conjugate singular points closest to 
the origin at z — ±27ri; therefore 



lim 



C,+l 



2tt 



(10) 



If we take into account that 

vr/2 



sin 



2j 6d6 



o 



VgXj + 1/2) 

2r(j + 1) 
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then we conclude that 



lim 




= lim 
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In other words, the p-power series has a finite radius of convergence R p = 2tx 
because of a pair of complex conjugate singular points at p c = ±2iri. 

Fig. [I] shows the first partial sums S? ] (p) = T + T lP + . . . + T N p N and 
the accurate numerical values of T(A). It dramatically illustrates the effect of 
the nonzero convergence radius Ra = \Z~2ir of the small-amplitude expansion 
determined by the singular points of T(A) closest to the origin in the complex 
^4-plane. 

We can provide another argument about the location of the singular points 
of T(A). First, note that the effective potential-energy function V(u) given 
in equation (J5J) exhibits a minimum V(0) = (1 — e p )/(2p) < and that the 
energy of the oscillatory motion is E = ii 2 /2 + V(u) = for the given initial 
conditions. Therefore, we expect a critical value of p given by V(0) = that 
yields p c = ±2iri in agreement with the analysis above based on the small- 
amplitude series. 

In order to verify those exact analytical results in a numerical way we con- 
structed Pade approximants [N, N] (p) [6] from the partial sums S^ N ^ (p) and 
looked for the complex zeroes of the denominator. A sequence of such zeroes 
appeared to converge to a limit quite close to ±6.3i with a small real part that 
was negligible compared to the errors of the estimates. Besides, assuming that 
there is an algebraic singular point [6] closest to origin of the form (z — zo) a 
we carried out the same Pade analysis, but now on T~ l dT / dA (as a function 
of p), and obtained roughly the same complex numbers that are quite close to 
±27ri. Therefore, there appears to be no doubt that the radius of convergence 
of the p-power series is in fact R p = 2n and is due to complex conjugate 
singular points located on the imaginary axis of the complex p-plane at ±27ri. 
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Recently, Amore and Fernandez [5] investigated the possible advantages of the 
inverted perturbation series that in the present case takes the form 



p = Pl AT + p 2 AT 2 + ... 

4AT AT 2 13AT 3 

~ ~ + 24vr 3 + " ' 

AT = T-2tt (13) 



Fig. [2] shows that the partial sums for the inverted series Sy*(AT) = p\AT + 
P2 AT 2 + . . . +pnAT n converge smoothly from below towards the accurate nu- 
merical values of T(A). We are presently unable to prove such most interesting 
feature of the inverted series rigorously. 

Summarizing: Earlier studies on the nonlinear oscillator ([1]) [2H3] have clearly 
shown that the period is finite for all values of the amplitude. However, they 
did not cast any light on the failure of the harmonic balance (with the ansatz 
x HB (t) = Acos(ojt)) that predicts a singularity for A = 2. In this paper 
we suggest that the harmonic balance may be reflecting the singular points 
that determine the radius of convergence of the small-amplitude series for the 
period. We have exactly calculated the location of those singular points and 
concluded that the series converge for < A < \/2tx. This result may help to 
understand similar difficulties in future applications of the harmonic balance. 
We expect that a harmonic-balance approach with more terms will give a 
frecuency with only complex singular points. 

In addition to what was mentioned above, we have shown that for this problem 
the inverted perturbation series appears to converge smoothly from below and 
it is therefore preferable to the original small-amplitude expansion. It is an 
old and well-known approach that may, in some cases, lead to surprisingly 
accurate results [5]. 
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Fig. 1. Accurate numerical period (solid line) and partial sums for the small-ampli- 
tude approximation (dashed lines). 



T 




Fig. 2. Accurate numerical period (solid line) and partial sums for the inverse per- 
turbation series (dashed lines). 



